Cryo-EM reveals binding of linoleic acid to SARS-CoV-2 spike glycoprotein, suggesting an antiviral treatment strategy.

The COVID-19 pandemic and concomitant lockdowns presented a global health challenge and triggered unprecedented research efforts to elucidate the molecular mechanisms and pathogenicity of SARS-CoV-2. The spike glycoprotein decorating the surface of SARS-CoV-2 virions is a prime target for vaccine development, antibody therapy and serology as it binds the host cell receptor and is central for viral cell entry. The electron cryo-microscopy structure of the spike protein revealed a hydrophobic pocket in the receptor-binding domain that is occupied by an essential fatty acid, linoleic acid (LA). The LA-bound spike protein adopts a non-infectious locked conformation which is more stable than the infectious form and shields important immunogenic epitopes. Here, the impact of LA binding on viral infectivity and replication, and the evolutionary conservation of the pocket in other highly pathogenic coronaviruses, including SARS-CoV-2 variants of concern (VOCs), are reviewed. The importance of LA metabolic products, the eicosanoids, in regulating the human immune response and inflammation is highlighted. Lipid and fatty-acid binding to a hydrophobic pocket in proteins on the virion surface appears to be a broader strategy employed by viruses, including picornaviruses and Zika virus. Ligand binding stabilizes their protein structure and assembly, and downregulates infectivity. In the case of rhinoviruses, this has been exploited to develop small-molecule antiviral drugs that bind to the hydrophobic pocket. The results suggest a COVID-19 antiviral treatment based on the LA-binding pocket.

Cryo-EM structure of a SARS-CoV-2 omicron spike protein ectodomain.

The omicron variant of SARS-CoV-2 has been spreading rapidly across the globe. The virus-surface spike protein plays a critical role in the cell entry and immune evasion of SARS-CoV-2. Here we determined the 3.0 Å cryo-EM structure of the omicron spike protein ectodomain. In contrast to the original strain of SARS-CoV-2 where the receptor-binding domain (RBD) of the spike protein takes a mixture of open ("standing up") and closed ("lying down") conformations, the omicron spike molecules are predominantly in the open conformation, with one upright RBD ready for receptor binding. The open conformation of the omicron spike is stabilized by enhanced inter-domain and inter-subunit packing, which involves new mutations in the omicron strain. Moreover, the omicron spike has undergone extensive mutations in RBD regions where known neutralizing antibodies target, allowing the omicron variant to escape immune surveillance aimed at the original viral strain. The stable open conformation of the omicron spike sheds light on the cell entry and immune evasion mechanisms of the omicron variant.

Conformational Flexibility and Local Frustration in the Functional States of the SARS-CoV-2 Spike B.1.1.7 and B.1.351 Variants: Mutation-Induced Allosteric Modulation Mechanism of Functional Dynamics and Protein Stability.

Structural and functional studies of the SARS-CoV-2 spike proteins have recently determined distinct functional states of the B.1.1.7 and B.1.351 spike variants, providing a molecular framework for understanding the mechanisms that link the effect of mutations with the enhanced virus infectivity and transmissibility. A detailed dynamic and energetic analysis of these variants was undertaken in the present work to quantify the effects of different mutations on functional conformational changes and stability of the SARS-CoV-2 spike protein. We employed the efficient and accurate coarse-grained (CG) simulations of multiple functional states of the D614G mutant, B.1.1.7 and B.1.351 spike variants to characterize conformational dynamics of the SARS-CoV-2 spike proteins and identify dynamic signatures of the functional regions that regulate transitions between the closed and open forms. By combining molecular simulations with full atomistic reconstruction of the trajectories and the ensemble-based mutational frustration analysis, we characterized how the intrinsic flexibility of specific spike regions can control functional conformational changes required for binding with the host-cell receptor. Using the residue-based mutational scanning of protein stability, we determined protein stability hotspots and identified potential energetic drivers favoring the receptor-accessible open spike states for the B.1.1.7 and B.1.351 spike variants. The results suggested that modulation of the energetic frustration at the inter-protomer interfaces can serve as a mechanism for allosteric couplings between mutational sites and the inter-protomer hinges of functional motions. The proposed mechanism of mutation-induced energetic frustration may result in greater adaptability and the emergence of multiple conformational states in the open form. This study suggested that SARS-CoV-2 B.1.1.7 and B.1.351 variants may leverage the intrinsic plasticity of functional regions in the spike protein for mutation-induced modulation of protein dynamics and allosteric regulation to control binding with the host cell receptor.

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex.

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

Receptor binding and complex structures of human ACE2 to spike RBD from omicron and delta SARS-CoV-2.

The coronavirus disease 2019 (COVID-19) pandemic continues worldwide with many variants arising, some of which are variants of concern (VOCs). A recent VOC, omicron (B.1.1.529), which obtains a large number of mutations in the receptor-binding domain (RBD) of the spike protein, has risen to intense scientific and public attention. Here, we studied the binding properties between the human receptor ACE2 (hACE2) and the VOC RBDs and resolved the crystal and cryoelectron microscopy structures of the omicron RBD-hACE2 complex as well as the crystal structure of the delta RBD-hACE2 complex. We found that, unlike alpha, beta, and gamma, omicron RBD binds to hACE2 at a similar affinity to that of the prototype RBD, which might be due to compensation of multiple mutations for both immune escape and transmissibility. The complex structures of omicron RBD-hACE2 and delta RBD-hACE2 reveal the structural basis of how RBD-specific mutations bind to hACE2.

E484K and N501Y SARS-CoV 2 spike mutants Increase ACE2 recognition but reduce affinity for neutralizing antibody.

SARS-CoV2 mutants B.1.1.7, B.1.351, and P.1 contain a key mutation N501Y. B.1.135 and P.1 lineages have another mutation, E484K. Here, we decode the effect of these two mutations on the host receptor, ACE2, and neutralizing antibody (B38) recognition. The N501Y RBD mutant binds to ACE2 with higher affinity due to improved π-π stacking and π-cation interactions. The higher binding affinity of the E484K mutant is caused due to the formation of additional hydrogen bond and salt-bridge interactions with ACE2. Both the mutants bind to the B38 antibody with reduced affinity due to the loss of several hydrogen-bonding interactions. The insights obtained from the study are crucial to interpret the increased transmissibility and reduced neutralization efficacy of rapidly emerging SARS-CoV2 VOCs.

Molecular insights into receptor binding of recent emerging SARS-CoV-2 variants.

Multiple SARS-CoV-2 variants of concern (VOCs) have been emerging and some have been linked to an increase in case numbers globally. However, there is yet a lack of understanding of the molecular basis for the interactions between the human ACE2 (hACE2) receptor and these VOCs. Here we examined several VOCs including Alpha, Beta, and Gamma, and demonstrate that five variants receptor-binding domain (RBD) increased binding affinity for hACE2, and four variants pseudoviruses increased entry into susceptible cells. Crystal structures of hACE2-RBD complexes help identify the key residues facilitating changes in hACE2 binding affinity. Additionally, soluble hACE2 protein efficiently prevent most of the variants pseudoviruses. Our findings provide important molecular information and may help the development of novel therapeutic and prophylactic agents targeting these emerging mutants.

Cryo-EM structure determination of small proteins by nanobody-binding scaffolds (Legobodies).

We describe a general method that allows structure determination of small proteins by single-particle cryo-electron microscopy (cryo-EM). The method is based on the availability of a target-binding nanobody, which is then rigidly attached to two scaffolds: 1) a Fab fragment of an antibody directed against the nanobody and 2) a nanobody-binding protein A fragment fused to maltose binding protein and Fab-binding domains. The overall ensemble of ∼120 kDa, called Legobody, does not perturb the nanobody-target interaction, is easily recognizable in EM images due to its unique shape, and facilitates particle alignment in cryo-EM image processing. The utility of the method is demonstrated for the KDEL receptor, a 23-kDa membrane protein, resulting in a map at 3.2-Šoverall resolution with density sufficient for de novo model building, and for the 22-kDa receptor-binding domain (RBD) of SARS-CoV-2 spike protein, resulting in a map at 3.6-Šresolution that allows analysis of the binding interface to the nanobody. The Legobody approach thus overcomes the current size limitations of cryo-EM analysis.

Potential Molecular Mechanisms of Rare Anti-Tumor Immune Response by SARS-CoV-2 in Isolated Cases of Lymphomas.

Recently, two cases of complete remission of classical Hodgkin lymphoma (cHL) and follicular lymphoma (FL) after SARS-CoV-2 infection were reported. However, the precise molecular mechanism of this rare event is yet to be understood. Here, we hypothesize a potential anti-tumor immune response of SARS-CoV-2 and based on a computational approach show that: (i) SARS-CoV-2 Spike-RBD may bind to the extracellular domains of CD15, CD27, CD45, and CD152 receptors of cHL or FL and may directly inhibit cell proliferation. (ii) Alternately, upon internalization after binding to these CD molecules, the SARS-CoV-2 membrane (M) protein and ORF3a may bind to gamma-tubulin complex component 3 (GCP3) at its tubulin gamma-1 chain (TUBG1) binding site. (iii) The M protein may also interact with TUBG1, blocking its binding to GCP3. (iv) Both the M and ORF3a proteins may render the GCP2-GCP3 lateral binding where the M protein possibly interacts with GCP2 at its GCP3 binding site and the ORF3a protein to GCP3 at its GCP2 interacting residues. (v) Interactions of the M and ORF3a proteins with these gamma-tubulin ring complex components potentially block the initial process of microtubule nucleation, leading to cell-cycle arrest and apoptosis. (vi) The Spike-RBD may also interact with and block PD-1 signaling similar to pembrolizumab and nivolumab- like monoclonal antibodies and may induce B-cell apoptosis and remission. (vii) Finally, the TRADD interacting "PVQLSY" motif of Epstein-Barr virus LMP-1, that is responsible for NF-kB mediated oncogenesis, potentially interacts with SARS-CoV-2 Mpro, NSP7, NSP10, and spike (S) proteins, and may inhibit the LMP-1 mediated cell proliferation. Taken together, our results suggest a possible therapeutic potential of SARS-CoV-2 in lymphoproliferative disorders.

A highly potent antibody effective against SARS-CoV-2 variants of concern.

Control of the ongoing SARS-CoV-2 pandemic is endangered by the emergence of viral variants with increased transmission efficiency, resistance to marketed therapeutic antibodies, and reduced sensitivity to vaccine-induced immunity. Here, we screen B cells from COVID-19 donors and identify P5C3, a highly potent and broadly neutralizing monoclonal antibody with picomolar neutralizing activity against all SARS-CoV-2 variants of concern (VOCs) identified to date. Structural characterization of P5C3 Fab in complex with the spike demonstrates a neutralizing activity defined by a large buried surface area, highly overlapping with the receptor-binding domain (RBD) surface necessary for ACE2 interaction. We further demonstrate that P5C3 shows complete prophylactic protection in the SARS-CoV-2-infected hamster challenge model. These results indicate that P5C3 opens exciting perspectives either as a prophylactic agent in immunocompromised individuals with poor response to vaccination or as combination therapy in SARS-CoV-2-infected individuals.

Cryo-EM analysis of the HCoV-229E spike glycoprotein reveals dynamic prefusion conformational changes.

Coronaviruses spike (S) glycoproteins mediate viral entry into host cells by binding to host receptors. However, how the S1 subunit undergoes conformational changes for receptor recognition has not been elucidated in Alphacoronavirus. Here, we report the cryo-EM structures of the HCoV-229E S trimer in prefusion state with two conformations. The activated conformation may pose the potential exposure of the S1-RBDs by decreasing of the interaction area between the S1-RBDs and the surrounding S1-NTDs and S1-RBDs compared to the closed conformation. Furthermore, structural comparison of our structures with the previously reported HCoV-229E S structure showed that the S trimers trended to open the S2 subunit from the closed conformation to open conformation, which could promote the transition from pre- to postfusion. Our results provide insights into the mechanisms involved in S glycoprotein-mediated Alphacoronavirus entry and have implications for vaccine and therapeutic antibody design.

Site-Specific Steric Control of SARS-CoV-2 Spike Glycosylation.

A central tenet in the design of vaccines is the display of native-like antigens in the elicitation of protective immunity. The abundance of N-linked glycans across the SARS-CoV-2 spike protein is a potential source of heterogeneity among the many different vaccine candidates under investigation. Here, we investigate the glycosylation of recombinant SARS-CoV-2 spike proteins from five different laboratories and compare them against S protein from infectious virus, cultured in Vero cells. We find patterns that are conserved across all samples, and this can be associated with site-specific stalling of glycan maturation that acts as a highly sensitive reporter of protein structure. Molecular dynamics simulations of a fully glycosylated spike support a model of steric restrictions that shape enzymatic processing of the glycans. These results suggest that recombinant spike-based SARS-CoV-2 immunogen glycosylation reproducibly recapitulates signatures of viral glycosylation.

Cryo-EM Structures of SARS-CoV-2 Spike without and with ACE2 Reveal a pH-Dependent Switch to Mediate Endosomal Positioning of Receptor-Binding Domains.

The SARS-CoV-2 spike employs mobile receptor-binding domains (RBDs) to engage the human ACE2 receptor and to facilitate virus entry, which can occur through low-pH-endosomal pathways. To understand how ACE2 binding and low pH affect spike conformation, we determined cryo-electron microscopy structures-at serological and endosomal pH-delineating spike recognition of up to three ACE2 molecules. RBDs freely adopted "up" conformations required for ACE2 interaction, primarily through RBD movement combined with smaller alterations in neighboring domains. In the absence of ACE2, single-RBD-up conformations dominated at pH 5.5, resolving into a solitary all-down conformation at lower pH. Notably, a pH-dependent refolding region (residues 824-858) at the spike-interdomain interface displayed dramatic structural rearrangements and mediated RBD positioning through coordinated movements of the entire trimer apex. These structures provide a foundation for understanding prefusion-spike mechanics governing endosomal entry; we suggest that the low pH all-down conformation potentially facilitates immune evasion from RBD-up binding antibody.

Structural Insight into the Binding of Cyanovirin-N with the Spike Glycoprotein, Mpro and PLpro of SARS-CoV-2: Protein-Protein Interactions, Dynamics Simulations and Free Energy Calculations.

The emergence of COVID-19 continues to pose severe threats to global public health. The pandemic has infected over 171 million people and claimed more than 3.5 million lives to date. We investigated the binding potential of antiviral cyanobacterial proteins including cyanovirin-N, scytovirin and phycocyanin with fundamental proteins involved in attachment and replication of SARS-CoV-2. Cyanovirin-N displayed the highest binding energy scores (-16.8 ± 0.02 kcal/mol, -12.3 ± 0.03 kcal/mol and -13.4 ± 0.02 kcal/mol, respectively) with the spike protein, the main protease (Mpro) and the papainlike protease (PLpro) of SARS-CoV-2. Cyanovirin-N was observed to interact with the crucial residues involved in the attachment of the human ACE2 receptor. Analysis of the binding affinities calculated employing the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) approach revealed that all forms of energy, except the polar solvation energy, favourably contributed to the interactions of cyanovirin-N with the viral proteins. With particular emphasis on cyanovirin-N, the current work presents evidence for the potential inhibition of SARS-CoV-2 by cyanobacterial proteins, and offers the opportunity for in vitro and in vivo experiments to deploy the cyanobacterial proteins as valuable therapeutics against COVID-19.

Possible Link between Higher Transmissibility of Alpha, Kappa and Delta Variants of SARS-CoV-2 and Increased Structural Stability of Its Spike Protein and hACE2 Affinity.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) outbreak in December 2019 has caused a global pandemic. The rapid mutation rate in the virus has created alarming situations worldwide and is being attributed to the false negativity in RT-PCR tests. It has also increased the chances of reinfection and immune escape. Recently various lineages namely, B.1.1.7 (Alpha), B.1.617.1 (Kappa), B.1.617.2 (Delta) and B.1.617.3 have caused rapid infection around the globe. To understand the biophysical perspective, we have performed molecular dynamic simulations of four different spikes (receptor binding domain)-hACE2 complexes, namely wildtype (WT), Alpha variant (N501Y spike mutant), Kappa (L452R, E484Q) and Delta (L452R, T478K), and compared their dynamics, binding energy and molecular interactions. Our results show that mutation has caused significant increase in the binding energy between the spike and hACE2 in Alpha and Kappa variants. In the case of Kappa and Delta variants, the mutations at L452R, T478K and E484Q increased the stability and intra-chain interactions in the spike protein, which may change the interaction ability of neutralizing antibodies to these spike variants. Further, we found that the Alpha variant had increased hydrogen interaction with Lys353 of hACE2 and more binding affinity in comparison to WT. The current study provides the biophysical basis for understanding the molecular mechanism and rationale behind the increase in the transmissivity and infectivity of the mutants compared to wild-type SARS-CoV-2.

Structures and distributions of SARS-CoV-2 spike proteins on intact virions.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virions are surrounded by a lipid bilayer from which spike (S) protein trimers protrude1. Heavily glycosylated S trimers bind to the angiotensin-converting enzyme 2 receptor and mediate entry of virions into target cells2-6. S exhibits extensive conformational flexibility: it modulates exposure of its receptor-binding site and subsequently undergoes complete structural rearrangement to drive fusion of viral and cellular membranes2,7,8. The structures and conformations of soluble, overexpressed, purified S proteins have been studied in detail using cryo-electron microscopy2,7,9-12, but the structure and distribution of S on the virion surface remain unknown. Here we applied cryo-electron microscopy and tomography to image intact SARS-CoV-2 virions and determine the high-resolution structure, conformational flexibility and distribution of S trimers in situ on the virion surface. These results reveal the conformations of S on the virion, and provide a basis from which to understand interactions between S and neutralizing antibodies during infection or vaccination.

Understanding structural malleability of the SARS-CoV-2 proteins and relation to the comorbidities.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a causative agent of the coronavirus disease (COVID-19), is a part of the $\beta $-Coronaviridae family. The virus contains five major protein classes viz., four structural proteins [nucleocapsid (N), membrane (M), envelop (E) and spike glycoprotein (S)] and replicase polyproteins (R), synthesized as two polyproteins (ORF1a and ORF1ab). Due to the severity of the pandemic, most of the SARS-CoV-2-related research are focused on finding therapeutic solutions. However, studies on the sequences and structure space throughout the evolutionary time frame of viral proteins are limited. Besides, the structural malleability of viral proteins can be directly or indirectly associated with the dysfunctionality of the host cell proteins. This dysfunctionality may lead to comorbidities during the infection and may continue at the post-infection stage. In this regard, we conduct the evolutionary sequence-structure analysis of the viral proteins to evaluate their malleability. Subsequently, intrinsic disorder propensities of these viral proteins have been studied to confirm that the short intrinsically disordered regions play an important role in enhancing the likelihood of the host proteins interacting with the viral proteins. These interactions may result in molecular dysfunctionality, finally leading to different diseases. Based on the host cell proteins, the diseases are divided in two distinct classes: (i) proteins, directly associated with the set of diseases while showing similar activities, and (ii) cytokine storm-mediated pro-inflammation (e.g. acute respiratory distress syndrome, malignancies) and neuroinflammation (e.g. neurodegenerative and neuropsychiatric diseases). Finally, the study unveils that males and postmenopausal females can be more vulnerable to SARS-CoV-2 infection due to the androgen-mediated protein transmembrane serine protease 2.

Mutation in a SARS-CoV-2 Haplotype from Sub-Antarctic Chile Reveals New Insights into the Spike's Dynamics.

The emergence of SARS-CoV-2 variants, as observed with the D614G spike protein mutant and, more recently, with B.1.1.7 (501Y.V1), B.1.351 (501Y.V2) and B.1.1.28.1 (P.1) lineages, represent a continuous threat and might lead to strains of higher infectivity and/or virulence. We report on the occurrence of a SARS-CoV-2 haplotype with nine mutations including D614G/T307I double-mutation of the spike. This variant expanded and completely replaced previous lineages within a short period in the subantarctic Magallanes Region, southern Chile. The rapid lineage shift was accompanied by a significant increase of cases, resulting in one of the highest incidence rates worldwide. Comparative coarse-grained molecular dynamic simulations indicated that T307I and D614G belong to a previously unrecognized dynamic domain, interfering with the mobility of the receptor binding domain of the spike. The T307I mutation showed a synergistic effect with the D614G. Continuous surveillance of new mutations and molecular analyses of such variations are important tools to understand the molecular mechanisms defining infectivity and virulence of current and future SARS-CoV-2 strains.

Structural basis for bivalent binding and inhibition of SARS-CoV-2 infection by human potent neutralizing antibodies.

Neutralizing monoclonal antibodies (nAbs) to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) represent promising candidates for clinical intervention against coronavirus disease 2019 (COVID-19). We isolated a large number of nAbs from SARS-CoV-2-infected individuals capable of disrupting proper interaction between the receptor binding domain (RBD) of the viral spike (S) protein and the receptor angiotensin converting enzyme 2 (ACE2). However, the structural basis for their potent neutralizing activity remains unclear. Here, we report cryo-EM structures of the ten most potent nAbs in their native full-length IgG-form or in both IgG-form and Fab-form bound to the trimeric S protein of SARS-CoV-2. The bivalent binding of the full-length IgG is found to associate with more RBDs in the "up" conformation than the monovalent binding of Fab, perhaps contributing to the enhanced neutralizing activity of IgG and triggering more shedding of the S1 subunit from the S protein. Comparison of a large number of nAbs identified common and unique structural features associated with their potent neutralizing activities. This work provides a structural basis for further understanding the mechanism of nAbs, especially through revealing the bivalent binding and its correlation with more potent neutralization and the shedding of S1 subunit.

Bat and pangolin coronavirus spike glycoprotein structures provide insights into SARS-CoV-2 evolution.

In recognizing the host cellular receptor and mediating fusion of virus and cell membranes, the spike (S) glycoprotein of coronaviruses is the most critical viral protein for cross-species transmission and infection. Here we determined the cryo-EM structures of the spikes from bat (RaTG13) and pangolin (PCoV_GX) coronaviruses, which are closely related to SARS-CoV-2. All three receptor-binding domains (RBDs) of these two spike trimers are in the "down" conformation, indicating they are more prone to adopt the receptor-binding inactive state. However, we found that the PCoV_GX, but not the RaTG13, spike is comparable to the SARS-CoV-2 spike in binding the human ACE2 receptor and supporting pseudovirus cell entry. We further identified critical residues in the RBD underlying different activities of the RaTG13 and PCoV_GX/SARS-CoV-2 spikes. These results collectively indicate that tight RBD-ACE2 binding and efficient RBD conformational sampling are required for the evolution of SARS-CoV-2 to gain highly efficient infection.